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The spin coherence phenomena and the possibihty of their observation in nanomagnetic insulators 
attract more and more attention in the last several years. Recently it has been shown that in these 
systems in large transverse magnetic field there can be a fairly narrow "coherence window" for 
phonon and nuclear spin-mediated decoherence. What kind of spin dynamics can then be expected 
in this window in a crystal of magnetic nanomolecules coupled to phonons, to nuclear spin bath 
and to each other via dipole-dipole interactions'? Studying multispin correlations, we determine the 
region of parameters where "coherent clusters" of collective spin excitations can appear. Although 
two particular systems, namely crystals of Feg-triazacyclonane and Mni2-acetate molecules, are 
used in this work to illustrate the results, here we are not trying to predict an existence of collective 
coherent dynamics in some particular system. Instead, we discuss the way how any crystalline 
system of dipole-dipole coupled nanomolecules can be analyzed to decide whether this system is 
suitable for attempts to observe coherent dynamics. The presented analysis can be useful in the 
search for magnetic systems showing the spin coherence phenomena. 

PACS numbers: 



X3 
O 

o 



> 
o 

oo 

o 

O 



I 

O 

o 



I. INTRODUCTION 

In the last decade the quantum tunneling phenomenon 
in nanomagnetic insulators has been attracting exten- 
sive interest. Many experiments have been done to 
study the tunneling relaxation in the ensembles of mag- 
netic molecules with central molecular spins S^^iSi^ 
These molecules couple to each other via dipole-dipole 
interactionSfSiSiL^ to phonons^iSiiSiiiii^ and to nuclear 
gpjj^g|5ji3ji4 rpj^g early study of these systems in a low 
temperature regime has been concentrated mainly on the 
incoherent tunneling in low transverse fields, when the 
magnitude of the ground state tunneling splitting 2Ao 
(produced by the tunneling between two potential wells 
separated by a barrier of magnetic anisotropy; Ao is the 
tunneling matrix element) is small in comparison with 
the parameters describing interactions with environment 
providing anomalously high decoherence. During last 
several years more attention has been paid to the spin 
coherence phenomenaii^*iSiiLi2iiS 

As it has been shown recently^S in nanomagnetic insu- 
lators in large transverse fields, where Ao{H^) increases, 
there can be a field region (" coherence window" ) in which 
both phonon and nuclear spin-mediated decoherence are 
drastically reduced (electronic decoherence in magnetic 
insulators is absent). The existence of such coherence 
window is important both for fundamental physics (at- 
tempts to find materials showing coherent spin tunneling 
phenomenon) and for quantum device engineering (at- 
tempts to make a solid-state qubit). 

At very low temperature each molecule with large cen- 
tral spin Si can be modelled as a two-level system (Ap- 
pendixIXIl. whose Hamiltonian Hi = — A^ff — ^iff oper- 
ates in a subspace of only two lowest states of Si. With 
Ai and ^i being the ground state tunneling matrix ele- 



ment and the longitudinal bias acting on i-th molecule, 
this two-state representation is valid only if A^ is small 
in comparison with the spin gap Eg to the next levels. 
For example, in two well-known central spin \S\ — 10 
systems, i^eg-triazacyclonane (Fes) and Afrii2-acetate 
(Mni2), this condition is met since Eg ^ 5 K and 
^ 11 K respectively while the values of a zero- field tun- 
neling splitting are ~ 10~^ K and ^ 10~^^ K. 

Suppose that molecules do not interact with each 
other. Then the central spin of any molecule can os- 
cillate between states | t) and | |) (IA3|I and this process 
is described by the probability Pxt(t) (|X4|) . If A >> ^, 
the amplitude of these oscillations is w 1. If the central 
spin of each molecule is also isolated from its nuclear sub- 
system and from the phonon thermostat, the tunneling 
oscillations, being coherent, can last for an infinitely long 
time. Interactions with the nuclear spin and the phonon 
thermostats lead to decoherence and, after the so-called 
decoherence time r^, coherence will be suppressed and 
oscillations will disappear. 

The decoherence " quality factor" , giving an estimation 
for the number of coherent oscillations in the system be- 
fore coherence will be suppressed, is ~ ^/"f4n where 
70 — h/{AoT^) is the dimensionless decoherence rate. 
The contributions to the decoherence time from in- 
teractions with the nuclear spins and phonons (r™ and 
T^^, respectively) arei^SiSi 
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— 3^coth(A„/fcBT), (1) 



where Eg is the half-width of the Gaussian distribution 
of the hypcrfinc bias energies; O^i is the Debye energy; 
and f2o ~ Eg is the energy of small oscillations in the 
potential wells. 
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The goal of the present work is to study the spin dy- 
namics in ensembles of dipole-dipole coupled magnetic 
molecules in the coherence window for the nuclear spin 
and phonon degrees of freedom at times t < Tnp, 

r„, = min{r;",rf}. (2) 

Namely, in this work we would like to study the inter- 
nal dynamics of a temperature equilibrated system, but 
not the dynamics induced by the artificial preparation of 
a system at t = 0, say, in state | tt ... t) (the anal- 
ysis of the latter problem will be presented separately). 
From now on, for the sake of brevity, the coherence win- 
dow will be called the NPC-window (nuclear spins and 
phonons coherence window). To illustrate the results, all 
particular calculations will be based on the parameters 
for two systems, namely, for crystals of -Fes a-nd Mni2 
molecules. 



II. HAMILTONIAN AND INTERACTIONS 

At very low temperatures a set of molecules with cen- 
tral molecular spins \Si\ = S coupled to each other via 
the dipole-dipole interaction can be described by the ef- 
fective Hamiltonian: 

H = ^(-A,ff - er^f ) + (3) 

i ij 

where and are the Pauli matrixes; is the tun- 
neling matrix element; and ^f" is the bias acting on i-th 
molecule from external and nuclear fields. The last term 
in Q describes the dipolar coupling between pairs of 
molecules, separated by distance \ fij \ = \ fi — rj\: 

where Ed = (/io/47r)g2^|5'2; ^o/47r = IQ-^N/A^ (in the 
SI system of units); ge is the electronic g-factor; and /is 
is the Bohr magneton. Note that Hamiltonian does 
not include the interactions with phonons and nuclear 
spins. Instead, the known results^fliSi for the phonon 
and nuclear spin decoherence rates, Eq.(^, will be used. 

Coherence window for the nuclear spin and phonon 
channels of decoherence opens up at high transverse 
fields, where the value of the tunneling splitting becomes 
large in comparison with the parameters describing in- 
teractions of the central spin with the environment. At 
these conditions all in a sample are approximately 
the same22 and for brevity can be replaced (where it is 
reasonable) by one parameter Aq, whose transverse field 
dependence Ao{H^) can be calculated using the corre- 
sponding molecular Hamiltonian for the central spin S. 
For both the Fes and the AIni2 molecules these Hamil- 
tonians are (approximately) known. 



(1) The "central spin" Hamiltonians for Fes 
and Mni2 molecules. Below ^ 10 K for Fes and 
below ^ 40 K for Mni2 these molecules are described by 
two similar 5* = 10 Hamiltonians of magnetic anisotropy: 

Hf"^ = -DSl + ESl + Ki{SX + St)~ get^sHS, (5) 

with24 D/kB = 0.23 K, E/kB = 0.094 K, and Ki/ks = 
-3.28 X 10-^ K; and 

^(Mn) ^ _r>sl-K^}st + Ki{Sl + St)-gefiBHS, (6) 

withes D/kB = 0.548 K, K^/kB = 1.17 x lO^^ K, 
Ki/kB = 2.2 X 10-^ K. 

Note that in the Fes system the tunneling splitting 
/S.o{F[^) and its period of oscillations with have 
been measured^SiSi while in the Mni2 system these pa- 
rameters have never been measured. The latter makes 
it rather problematic to verify the value of the tunnel- 
ing splitting obtained directly from the Hamiltonian 10 . 
However, we would like to study the region of large trans- 
verse fields where Ao(-ff^) is already large (although 
Ao << Eg) and is less sensitive to some variations of 
the anisotropy constants^S. (moreover, at some stage we 
start to make estimations rather than exact calculations) . 
Thus, in what follows we use the Hamiltonians Q and 
© for Fes and Mni2 molecules. 

(2) Dipolar interactions. For the sake of defi- 
niteness we apply a transverse magnetic field along the 
X-axis, so that only the and the Sf projections of 
the total molecular spin Si are nonzero. Therefore, the 
interaction term Vdd can be rewritten as: 

Vd4n,)^ E ^;/(^L)^r^/, (7) 

{a,P} = {x,z} 

where all V^f{fij) can be obtained from Eq.Q. The 
i-th bias energy ^f" in Q, as it is written, contains con- 
tributions only from the longitudinal external and nu- 
clear fields. The dipolar contribution to the total bias 
acting on z-th molecule can be written in the form 
= —gelJ.BSiH?{dip) and the longitudinal dipolar filed 
Hf{dip) at i-th site is: 

where Ed — (/io/47r)ge/iB'5' and z^ , Xij are the corre- 
sponding components of vector . 

The distributions of the dipolar bias energies created 
by molecular spins in polarized and depolarized samples 
are different in the low transverse field limit and similar 
in the high transverse fields limit (where Si is oriented 
nearly along the transverse field direction). At low trans- 
verse fields the half-width Wd of the dipolar bias dis- 
tribution in a completely depolarized sample is several 
times larger than in a polarized sample. At high trans- 
verse fields the half-width Wd in both samples is nearly 
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the same (comparing the longitudinal field distributions 
for polarized and depolarized samples at — 4.8 T in 
FigElof Appendix IbI one can see that they are nearly the 
same) . This parameter can be calculated numerically for 
any sample. 

(3) Hyperfine interactions. The interactions 
of the central molecular spin Si with the nuclear spin 
bath lead to the "spread" of each molecular spin state 
characterized by the half-width Eg of the Gaussian dis- 
tribution of the hyperfine bias energies ^at . For Nn nu- 
clear spins Ik in each molecule, one findsiSSiSiiSi E'^ = 

Y^k=ii^k + l)-^/£ (^11)^/3, where {lu^^} are the (longitudi- 
nal) couplings between the central spin and each k-th 
nuclear spinii^ii4 Knowledge of all nuclear moments and 
positions of all nuclei in the molecule^ allows one to cal- 
culate all these coupling constants and Eq^^^'^^'^^^ 

(4) The transverse magnetic field behavior of 
important parameters. The ground state | f|~) 
(symmetric) and the excited state | JJ-) (antisymmet- 
ric) ljA2p of Hamiltonian (|A1|) are separated by the en- 
ergy gap 2ei = 2(Af -I- Cf)^^'^- At low temperatures 
{keT < Ao) in the limit >> S,i in a temperature 
equilibrated sample most of molecules are in states | ff). 
Then, calculating matrix elements (fl^ |f/| = ^i/si and 
iil Iff I i[) = AJe,, for A, » one finds \Sf\ i[) ^ 
and (-|| \Sf\ ii) — > S. Thus, as Ao increases with the 
transverse field, both Wd and Eg should decrease. 
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FIG. 1: Fes crystal (cluster of 50^ Fes molecules, for de- 
tails see Appendix 0. The curves are: circles - half- width 
Wd of the dipolar bias distribution vs transverse field (in 
Tesla) for model depolarized sample (see the text); stars - 
Wd{H^) for polarized sample; triangles - same as circles, but 
for the field-compensated sample; diamonds - half-width Eo 
of the hyperfine bias distribution; squares - same as diamonds 
but for the field-compensated sample; solid line - tunneling 
matrix element Ao{H^). All molecules are in a state | if). 
(Wd, Eo, Ao are in Kelvins.) 

When Ao << Wd, in a depolarized sample the value 
of Wd is several times larger than in a polarized sample. 



To understand how Wd behaves at large transverse fields, 
it is sufficient to calculate this parameter in a model de- 
polarized sample where all molecules are in states | ff), 
but Sf /\S^\ = Om^ The transverse field dependence 
of important parameters for crystals of Feg and Mni2 
molecules is presented in Figs^ and [3 (the description 
of our calculation procedure is given in Appendix IB|) 
Deviations from the results of Figs^ and [3 for nonzero 
populations of states | -IJ.) are insignificant up to the limit 
of equipopulation - this is clear, for example, from Fig[Sl 
(Appendix • 

Depending on the crystal structure and the sample 
geometry, the dipolar fields distribution can be shifted 
(such a shift can be rather large, see, for example, Fig[Sl 
in Appendix IB|) . This shift changes with the transverse 
field. The larger the shift, the slower both the Wd{H^) 
and the Eo{H^) decrease with H^. This effect can be 
seen in the high-field part of Fig^ Two upper curves 
for both Wd and Eo represent the results of calculations 
in our Feg cluster "as it is" (with no longitudinal field 
compensation). To obtain the two lower curves for both 
Wd and Eg, the corresponding external longitudinal field 
was applied for each value of the external transverse field 
to shift a position of the longitudinal fields distribution 
back to zero (the longitudinal field compensated sample). 
The shift in our Mni2 sample is small. 

1x1 oS (Mn-i 2 crystal) 




FIG. 2: Mni2 crystal (cluster of 50^ Mni2 molecules; no 
"faster relaxing species"—*—). The curves are: solid line with 
circles - WoiH'^) for "depolarized" sample; stars - WoiH"^) 
for polarized sample; diamonds - Eo{H^); solid line - Ao{H^). 



(5) NPC-window. When studying the spin dynam- 
ics in the NPC-window, one needs to know the region of 
the field where this window is situated. This window can 
be rather narrow and for our examples of the Fes and 
the Mni2 systems this can be seen in Fig|31 Since we 
are not going to discuss here a coherence optimization 
strategyjSS in this Figure we present the transverse field 
behavior of the dimensionless decoherence rates 7^" and 
Eq.J^I, for external field along the x axis only and 
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for molecules containing only natural isotopesi^ (Note 
that both rates are almost insensitive to changes in the 
populations of states | ff) and | -IJ-).) The small oscillation 
energy fio is ~ Eg and, like Eg{H-^), slowly decreases 
with H-^. In zero field Qo = 2SC^{DEf/'^ (C_l « 1.56) 
for Feg and Vlo ^ 2SD for Mni2^ The Debye energy 
Qjj for Fes ^^nd Mni2 is known experimentally^ 
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FIG. 3: The transverse magnetic field behavior of the di- 
mensionless decoherence rates 7^" (shown by squares and 
circles) and 7^'' {diamonds and triangles) at azimuthal an- 
gle ifi = (i.e., along the x-axis) for Fes {squares and di- 
amonds) and Mni2 {circles and triangles) systems at T = 
0.05 K. The results are presented for the natural isotopes 
{Fe^^, H\ Br''^, iV", C^^ and O^*^ species for Feg molecule 
and Mr^5^ H^, C^^ and O^^ species for Mni2 molecule). 



III. MULTIMOLECULAR PROCESSES IN THE 
LIMIT Ao » {WD,Eo,h/T„p} 

If in the NPC-window the half-width Wd of the dipo- 
lar bias distribution is larger than Aq, any spin dynamics 
is, in general, incoherent. Outside of this window, inde- 
pendently of the ratio Aq/Wd, the spin dynamics is also 
incoherent. In this Section we study the niultimolecu- 
lar correlations induced by the dipole-dipolc interactions 
between molecules in states | ff) and | JJ.) in the NPC- 
window. We assume that in the region of fields of our 
interest Aq >> {Wd, Eo,h/Tnp} (like, for example, in 
the region 3.5 < < 4.2 T for Fes and in the region 
6.7 < 7?^ < 7.4 T for Mni2, see FigEJ. 



A. One-pair processes 

At very low temperatures, when only two lowest states 
of each molecule are occupied, in the limit Ao >> 
{Wd, Eojh/Tnp} we work in the representation (|A2|) 
{| tr) J -U-)} of the Hamiltonian (|Xl1) ( Appendix EJ. Con- 
sider one pair of interacting molecules in the sample (Ap- 



pendix [nj. Such a pair is described by the Hamilto- 
nian (|C3I) {£,i is the bias energy; in general, it is the 
time-dependent parameter) and can be found in four 
states: | fTifTa); | ilih); I \ ^2)- In the limit 

Ao >> \V^f\ two central "flip-flop" states | friJJ.2) and 
I 4J.i'fl"2), linked by the effective tunneling matrix element 



A// - \V,^,%R)\ 



A1A2 

£l£2 



(9) 



(only the largest term in B{R) Ea. ljC6|) is shown - its 
form is similar to that known from the theory of dielectric 
glasses^&2i), can be considered as an effective two-level 
system with the asymmetry 

?// « |£i - £2! ^ le? - C2I/2A0 « Ao (10) 

(recall that all A^ are supposed to be the same). Note 
that all V^f {R) are independent of the external field. 

Two other states are separated from the two flip-flop 
states by the energy gaps > Aq and in the region of fields 
where Ao > jV"^"^^!, their effect on the flip-fiop transitions 
is small. If Afj > i^yj, two flip-flop states are in reso- 
nance and the amplitude of oscillations (with frequency 
Eff (?// + A/^)^/^ (EUl) between them is ?» 1. At 
the same time, transitions between other states are not 
in resonance (they are accompanied by the energy change 
> Ao) and their amplitude is < {V^f )'^ / Af . 

In the field region where Aq ^ I^Jd'l it is more con- 
venient to solve the problem for the dynamics of a pair 
of interacting spins in the basis set (|A3p .'^^ However, in 
this limit it can be rather difficult (or even impossible) 
to observe coherent dynamics in an ensemble of spins. 
First of all, in this limit Wd also can be > Aq (like in 
Fes and AIni2). Moreover, the variety of different collec- 
tive processes leads to additional phase randomness and, 
consequently, to suppression of coherence. 

In what follows we suppose to work only in the part of 
the NPC-window where Ao > \V2f\ and the probability 
to observe coherent spin dynamics is larger. 

At low temperatures {ksT < Ao) a number of the 
molecules in the excited state | JJ.) with energy -\-e can be 
estimated as Nex{T) ~ 7Voe-^/''«^/(e-'^/'=«^-he^/'=^'^) ~ 
A'oe"^'^"/'^^-^ {No is the total number of molecules) and 
is small compared to the number of the ground state 
molecules. These excited molecules are uniformly dis- 
tributed over the sample, and each of them is surrounded 
by the No/N^x ground state molecules. The ex- 
cited molecule can be in a flip-flop resonanc^^ with 
the ground state molecule only if C// ^5 ^ff- The 
time needed for the flip-flop transition to happen is 
~ fi/Aff ^ 0{R^) and the fastest transitions are ex- 
pected to be between the nearest-neighbor molecules. 
For the effective two-level systems composed of the two 
nearest-neighbor molecules we introduce the correspond- 
ing effective tunneling matrix element A^j and the asym- 
metry SJj^. 

In a simple cubic lattice each excited molecule can 
make a flip-flop transition with any of its six nearest- 
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neighbor ground state molecules with the same probabil- 
ities. In a generic lattice these probabilities can be dif- 
ferent since Ajj depends on the lattice structure. The 
average over three crystallographic axes value of the AJJ 
is ^ WoiH-^ = 0) for polarized sample, see Figs^-Ill 
{Wd{0) for polarized sample is ~ ED/Vj^\ where V^' 
is the volume per one molecule). 

If AJj > ^jjf for the overwhelming majority of the 
nearest-neighbor molecules (this issue is discussed in Sec- 
tion it is unlikely that at low temperatures any 
resonant pair of the nearest-neighbor molecules (say, i- 
th and j-th molecules) will remain in resonance for a 
long time. Instead, since the total probability for the 
excited molecule (either i-th, or j-th, as a result of oscil- 
lations ljC7|) ') to create a resonance with one of the other 
five nearest-neighbor molecules is larger than the prob- 
ability to remain in resonance with the same molecule 
all the time t ~ r„p, the fastest flip-flop transitions can 
"propagate" through the crystal involving more and more 
new molecules. Of course, not only the nearest-neighbor 
molecules can be involved, but also the "lengthy" pairs 
(with Aff{R) < Ayj). However, flip-flop transitions be- 
tween the nearest-neighbor molecules are faster. 

In what follows, for brevity, these "mobile" (or "poten- 
tially mobile" ) flip-flop transitions between the states | f|~) 
and I J|) in the nearest-neighbor molecules will be called 
flipons^^ (a kind of magnon). The number of flipons is 
determined by the number of excited molecules N^xiT). 
In a generic lattice A'^j can be different along different 
crystallographic axes. However, if flipon moves along one 
axis and for this axis {(A^j)^} are large in comparison 
with {(^^™)j}, such a movement is, in some sense, "co- 
herent" since flipon leaves site i only because of equal 
probabilities for the excited spin to create a resonance 
with both of its nearest neighbors along this axis. 

It is worth mentioning that, if there is a whole distribu- 
tion of Ao (say, if there are "faster relaxing species " ^^i^^ ) 
in a sample, the fraction of resonant flip-flop molecules 
decreases. This is because in such a sample for some 
fraction of pairs the asymmetry £^ff can be ^ Aq (and 
Ao increases with H^). These impurities can essentially 
limit (or even completely block) the motion of flipons. 

B. Multi-pair processes 

On average, two nearest-neighbor excited molecules are 
separated by the distance 

i?e.(T) ^ ivy^NjN,xy/\ (11) 

At kgT « Ao, Rf.x{T) is large compared to a = 

{V^')^^^ (in a cubic lattice a = a; a is the lattice 
constant). Consider two pairs of the resonant nearest- 
neighbor molecules and let the distance between these 
pairs be R. This group of four molecules contains two 
excited molecules (with energies £2, £2) ^^'^ two ground 
state molecules (with energies £1, Si). In the limit 



Ao >> {Wd,Eo, h/Tnp} and Ao > \V^f\ both these res- 
onant pairs experience mainly the flip-flop transitions. 

If i? >> a, the strength |V'^'^'^(i?)| of the interac- 
tions between molecules belonging to different pairs is 
< Ay J. Then, using the same arguments as for one pair 
of molecules, in the case of two resonant pairs we can also 
consider only corresponding collective "flip-flop" transi- 
tions between the eigenstates of each resonant pair. The 
effective tunneling matrix elements connecting these col- 
lective flip-flop states (separated from all other states by 
the energy gaps > AJj) is 

AfjiR)^\V,^,^R)\T&j:&, (12) 
^ff ^ff 

where E]f - iiCff? + {AJ^ff/^. Similarly to the case 
of one pair of molecules, these collective flip-flop states 
can also be considered as an effective two-level system 
with the asymmetry 

Note that we deliberately consider two pairs of nearest- 
neighbor molecules since transitions between such 
molecules are faster, and in the limit of our interest the 
probability to flnd them in resonance is larger. 

The effective matrix element A^^ describes the flip- 
flop transitions between the eigenstates of two resonant 
pairs (i.e., of two effective TLS). If A^^ > ^/J, two reso- 
nant pairs are in resonance with each other. The flip-flop 
transitions between the eigenstates | ff) and | JJ.) of each 

molecule inside of one resonant pair are described by the 

(2) 

matrix element A^J. Then, since A]] > Ayi, the fre- 
quency of oscillations between states | ff) and | ^) of 
each molecule in such a resonant group of four molecules 
is ^ ^//i but the group correlation time is ^ h/Aj^. 
In a generic lattice, if the nearest-neighbor molecules in 

two pairs are located along different axes, the asymmetry 

(2) (2) 

can be ~ A^J » A^^ and such two pairs can be 

out of resonance. However, if molecules in both pairs are 

located along the same axis (with a lattice constant a) 

and if AJ'j >> the asymmetry is 

4/ - I^Z-f («)l/(^^/Ao)' m/Ao)-Oi^f/Al), (14) 

where the average value of /(^i/Ao) can be estimated 
roughly as ~ W^/A^. Then, for the average asym- 

metry one gets ^/j ^ {V^i {a)\{WD / A^f « \V^^i{a)\ 

(2) (2^ 

and for such pairs the condition Ay^- <^ can be, in 
principle, fulflUcd (actually, the difference of two mean 
energies A and A (jC6|) also contributes to 4^'' and 
this gives a similar effect). The term neglected in (|14f) 
is ~ (AJ'J/Wd(H^))'^ times smaller than the retained 
one - in the field region of our interest in most systems 
A]]{a)/WD{H^) » 1. 

Note that for resonant pairs composed of the flip-flop 
molecules with A//(i?) << A^,, the glass-like scenario^ 
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can be realized. In this case two pairs can be in resonance 

(2) (2) 

only if Ajj ~ A// (for most of such pairs ~ A//). 
Knowing the sample average value of the asymmetry 

— (2) "* (2^ 

, from the requirement A\^{R) > one can esti- 
mate the average "resonant" distance between two pairs: 

Rlll{H^)-[V}'^V,,/^fj{H^)V^\ (15) 

where^? Vdd ^ Ed/V^'^\ If i? < R?X, two pairs could 
be, in principle, in resonance with each other. However, 
even if 0, not any two pairs are in resonance since 

if i? > Rph, where 

R,h{T,H^) ^ [v!;^^VMrf{T,H^)/h]"\ (16) 

the two-pairs correlation time is longer than incoherent 
phonon-assisted transitions in each molecule. Only the 
pairs satisfying the condition R < R,n = min{i?res, Rph}, 
can be in resonance. Thus, if Ri,j.{T, H^) < Rm, most of 
the closest pairs of resonant molecules are able to come 
into resonance with each other. This happens at temper- 
atures T > Tm, 

Tm --max{T^f„Tph}, (17) 

(2.) 

with Tres and Tph given by {ksT < A^): 

ksT^HiH^) ^ 2Ao/ \n[Vdd/^'j], (18) 

kBTph{H^) - 2AJ \n[Q%Vdd/S''nlAl]. (19) 
If at these conditions T„p > tc, where 

^ h/\vii{R,,)\ ^ in/Vdd){RlJvP), (20) 

at t > the whole hierarchy of (more or less) correlated 
flip-flop clusters of the increasing "size" n (the number 
of involved resonant flip-flop pairs) can, in principle, ap- 
pear. The time tc estimates the cluster correlation time. 

Note, however, that if A^™ > for most of the 
nearest-neighbor molecules, instead of interactions be- 
tween fixed pairs of resonant molecules, in the limit of 
our interest we have a set of flipons moving through the 
sample and interacting with each-other. At T < Tm they 
participate in the collective processes very rarely (interac- 
tions can be neglected). When temperature increases, the 
number of flipons also increases and collective processes 
become more frequent. At T > Tm correlations between 
flipons, in principle, may still lead to the creation of cor- 
related clusters. However, due to various decorrelation 
(dephasing) processes, these clusters can be destroyed 
rapidly (or they will not be able to appear at all). Such 
dephasing processes will be considered in Section Fill CI 

C. Decorrelation 
(1) Flipon motion. If flipons are delocalized, the 

(2) 

effective tunneling matrix element A , / changes with the 



distance between flipons resulting in the suppression of 
correlations in clusters (if they appear). 

Suppose that at i = there is a correlated cluster 
{Tm < T < Ao/ks) composed of nearly equidistant 
flipons (with distance w Rex)- If A^™ >> and if the 
flipons move along the same axis, correlations between 
them will not necessarily be destroyed immediately af- 
ter the flrst "jump" . Of course, if at i > the flipons 
start to move along different axes, in a generic lattice 
any correlations can be destroyed almost immediately 
(i.e., at < ^ h/A^"^) since in this case ^j-y can become 

(2) 

A™j >> A\j. Note, however, that the flipons have 
larger probability to move along the axis with shortest 
lattice constant. Then, if we consider only a quasi-ld mo- 
tion of flipons (i.e., along the same axis), we can estimate 
the longest " motional" dephasing time t™ . Comparison 
of this time with the cluster correlation time tc {tc < Tnp) 
shows whether the correlated cluster with the average 
distances Rex between the nearest-neighbor flipons can 
appear. 

For the sake of simplicity, we approximate the flipon 
centers of mass motion by the discrete "random walks" 
model (Appendix 0. At t > the distances R{t) be- 
tween the nearest-neighbor flipons in the whole cluster 
become distributed around Rex with nonzero half-width 

f 2) 

SR{t). Thus, instead of a single "line" Ajl{Rex) one also 

f 21 

gets a whole distribution of values A^-j (i?) with nonzero 

mean-square deviation 5A^^^\t) = {{Af})^ - (Ag)2)i/2. 

Knowing (5A(^-'(i), the motional dephasing time can be 
obtained from the condition 

NT 

rT=tfNT; 5]t/JA(^)(^)~n; tf^h/A]] (21) 

for N = t/tf (or from the condition /J"" dt 6A'^^'>{t) ~ h 
at large values of N) . Obviously, the correlations in the 
whole cluster will be destroyed together with the destruc- 
tion of resonances between the nearest-neighbor flipons. 

f 21 

Since in each pair both flipons can move, for {Aj^) we 
have 

N 

- E \Vii{Rex + {ri^r2)a)\PN{ri)PN{r2) (22) 

ri ,r2 — — N 

with the condition {Rex + (^i — r2)d) / Rex > vi^)- Here 
77 (T) is the dimensionless (in units of Rex{T)) minimally 
possible distance between the centers of mass of two 
flipons. Each distribution PN{ri) gives the probability 
to flnd the z-th flipon at the distance rid {ri < N) from 
its t = position after total N steps f Appendix ID|| .'^^ 

The solution of Eq.||2U depends on 77 and p (Ea. l\D4\i ) 
and can be found numerically. For p = q = s = l/3 
(Eq.(inill) and p = 2/3 we get 

rT = N^Hf EE 3A(r,)[i?e./S]'i/ - 3X{fj)[a/ Rex]tc. (23) 
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If ?7 = 2a/Rex, for 0.05 < 77 < 2/3 we get A = A2 « 
1.6r;2 + 0.3577 + 0.045. If = a/Re^, for 0.05 < < 
1/2 we get A = Ai « 0.3^ + 0.2?/ + 0.03. Note that 
the configurations of the nearly equidistant flipons with 
Rex = a do not exist, in contrast to those with R^x = 2a. 
However, if flipons move along the same axis, but in the 
nearest- neighbor rows, the centers of mass of some flipons 
can be separated by the distance a. To take this effect 
into account, one can use A = (A2 + Ai)/2 for estimations. 

The answer for can be found in the equivalent di- 
mensionless form DfT^/R^^ = X{ri,p), where Df is the 
flipon effective diffusion coefficient (|D4|I and at p = 2/3, 
A(77, p) = A(77) from Ea. l|^ . The p-dependence of t™ 
is roughly ^ p~^l^ . Then, t"^ jtc can become larger ei- 
ther (i) at r — > lS.o|kB^ when flipons are in their dense 
phase and essentially localized {SR{t) 0); or (ii) if 
s = l— p — g— >l,p— >0 and flipons are almost immo- 
bile even at T < Ao/ks- The latter can be, in principle, 
realized in a sample with impurities. 

If tJP{T) « tc{T), the creation of a correlated cluster 
at the average distance Rex{T) is virtually impossible. 
Solving either the equation N™{Ri,)tf = tc{Rc), or the 
equation 

i?/ie(i?c)/i?e=A(77e,p), (24) 

one finds the average distance Rc and the temperature 

ksT, ^ 2AoiH^)/ HiRc/af - 1] (25) 

at which t™ ~ tc and cluster can appear. For p — 2/3 
and = 2a/ Rc we get Rc ^ 3a and fc^Tc ^ 0.6Ao(-ff^). 
For p — 2/3 and 77c — a/Rc we get Rc ^ 2a and 
Tc ^ Ao{H^). These estimations shows that, if the sce- 
nario with p«g«s«l/3is realized, r™(T) < tc{T) 
almost everywhere except the flipons dense phase at 
T Ao/ks, where t^{T) - tc{T) and where tdT) de- 
creases itself (as well as Tnp(T), see Eq^. If, in contrast, 
s ^ 1 and p — > 0, t™ » tc and correlated clusters can 
appear even at Ts/ < T < Tc- 

(2) "Spectral diffusion". The above described 
picture is valid only if A > ^ for most of the nearest- 
neighbor molecules. In the opposite limit, at T > Tm 
the correlated clusters will be composed of almost immo- 
bile "lengthy" flip-flop pairs with A//(i? > a) < A^^, 

satisfying the condition A/y ~ ^// - This scenario is 
very similar to that in dielectric glasses^ and the clus- 
ter dephasing time at i < T„p will be determined by the 
process similar to the "spectral diffusion" in glassesi^^i^ 
The change of states of fixed effective TLS results in the 
bias fluctuations and, consequently, in the dephasing. In 
this hmit {AJ^ < ^^y) the cluster "spectral diffusion" 

(2) 

dephasing time rj is ~ tc since the asymmetry for 
most groups of two resonant pairs is ~ A/y . 

In the limit A^J > the spectral diffusion-like pro- 
cess can contribute as well. Instead of going deeply 
into the details, here we only estimate the correspond- 
ing effects. Consider, for simplicity, the case of immo- 
bile flipons. The bias (^f ? )», acting on any i-th flipon. 



contains contributions from all individual spins in the 
sample. When any j-th flipon makes a transition, the 
change of the bias, acting on ?-th flipon, is 5{j{Rij) ~ 
VII{R^J) {i]fm]f)j/{E]fUE]f),. (Here {^Jf).,, can 
be both positive and negative and the term VJJt^Tj does 
not change its sign when j-th flipon makes transition.) 
Then, if Nt flipons make a transition (max(7Vt) = Ngx), 
the total change of the bias acting on i-th flipons is 

Depending on the degree of "polarization" Mt = 
{Ng, - Ncs)/{Ngs + Ncs) of the group of Nt flipons 
(Ngs and Ncs are the numbers of flipons in their ground 
and excited states), the total bias change |i5/| can vary 
roughly from ?/^~ \ViiiRcx)\ [{^]fh/ {A]jy,][^]pA]J] 
to - S{ In(iVt) [A]] is the sample average value of ( AJJ)j 

and is the sample average absolute value of (^^j)j). 
Then the shortest dephasing time rj is ^ tf. However, 
the dipole-dipole interaction changes its sign with the 
direction of Rij and, on average, in the case Mt 1 
the "surface" spins will mainly determine the maximum 
value of the bias change \S(\. For spins (flipons) in the 
bulk this essentially reduces \6{\ and increases rj. 

If Mt 0, simultaneous transitions of many flipons 
nearly cancel the effect of each other resulting in |i5/| < 

S{ « |y//(i?e£c)| << Aj^. In this case the spectral dif- 
fusion mechanism cannot destroy the resonance neither 
inside of the individual flipons, nor between them. In- 
deed, the contribution from the to the asymmetry 

^fj in the limit AJ] > is ~ iCfff/2A]] and its 

change due to S{ is ~ S{ (^}']?)V(A^p^ 

On average, in a temperature equilibrated sample, it is 
plausible to assume Mt and in the limit A^™ > 

the spectral diffusion effect is much weaker than the mo- 

(2) 

tional dephasing effect (broadening of the Aj:j distribu- 
tion due to the flipons motion). This remains valid also 
if the flipons are allowed to move. In this case the spec- 
tral diffusion dephasing time rj is roughly ~ [AYj /Cf^Y 
times longer than the motional dephasing time r™ (as it 
was already noted, if flipons move along different axis, in 
a generic lattice correlations can be destroyed already at 
t^tf). 

IV. DISCUSSION 

In this Section we discuss a temperature equilibrated 
sample at ksT Aq, where only the fraction Nex(T) of 
molecules are in the excited states | ^). Thus, we are 
not going to discuss the problem of the magnetization 
relaxation in the limit Ag >> {Wd, Eq}- If, for example, 
the sample is prepared at very high temperatures and 
then rapidly cooled down to low temperatures, it will 
start to relax to its temperature equilibrated state - this 
relaxation process will not be discussed here either. 

If, in some system, in all the NFC-window Aq ^ IKzrf'l' 
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it can be very difficult (most probably, impossible) to 
observe any coherent spin dynamics in an ensemble of 
interacting spins f Section IIII A|l . Here we consider only 
the part of the NPC-window, where Aq > |V"j^^|. In both 
the -Feg and the Mni2 systems, which are used in this 
work to illustrate how the problem can be analyzed, this 
field region is situated to the right of the minimum of 
7™ + 7^^* (see FigOJ. The collective processes, having 
the largest amplitude in this region of fields, are the flip- 
flop processes. 



A. The ratio A]]/^ 



nn 

fS 



The fastest pair flip-flop processes are the processes 
between the nearest-neighbor molecules. The average 
strength of the nearest-neighbor dipole-dipole interac- 
tions and the average value of the flipon effective tun- 
neling matrix element Aj^ are ^ Vdd {Vdd ^ 0.12 K 
for Feg and ~ 0.07 K for Mriia). To find the aver- 
age asymmetry ^^y, we first calculate the distributions 

Pi2ie,H^) of the e = ei- 62 (e. = (A| + ^fy/^; 
all Ai depend on both the external and the dipolar 
transverse fields and are obtained by the exact diago- 
nalization of jSJ and (|SJ), Appendix IB|| . where £1^2 are 
for the nearest-neighbor molecules only. Then, we get 
= (1/2) J^^° de \e\Pi2{£,H^) {so is the maximum 
value of e). 
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nearest-neighbor molecules 



5000 



(b) Mni 2 crystal 
Hx=7.0 T 
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FIG. 4; The distributions of the e = e\ ~ 62 for the nearest- 
neighbor molecules: (a) at — 3.8 T in the crystal of Fes 
molecules (no field compensation); (b) at — 7.0 T in the 
crystal of Mni2 molecules (no faster relaxing species). All 
cluster parameters are the same as in Figs0-|3 Here e is in 
Kelvins. The states | ff) and | i}-) are equipopulated. Both 
distributions are averaged over three crystallographic axes. 

For example, two distributions (at — 3.8 T for Fes 
and at = 7.0 T for Mni2) are presented in Fig^ For 
these distributions we get « 1.5 x 10~^ K for Fes 
(with no field compensation) and £^Yf ~ 2.3 x 10^"* K for 
Mni2 (the corresponding mean-square deviations 5e = 
(e^ - are w 4.7 x lO'^ K and « 6.4 x 10"* K 



respectively). For the field-compensated Fe^ sample, the 
corresponding values are about two times smaller. 

To summarize the results, we note that under the as- 
sumption of the absence of impurities with larger (or 
smaller) values of A^, in both systems in the NPC- 
window the average asymmetry ^j?^ (as well as 5e) is 

small in comparison with Ayj for most of the nearest- 
neighbor molecules. In this case fiipons can move in both 
systems. 

Note, however, that in the Mni2 crystal there are 
faster relaxing (minor) specie3^2i2^ (about 5 — 10 % of 
all molecules have lower potential barrier and larger val- 
ues of Aj) - these species can change somehow both the 
and the Se. In the NPC-window for major species 
the spin dynamics of these minor species is incoherent al- 
ready. Since even the Hamiltonian for these minor species 
is still unknown (or unpublished), we cannot describe 
their effect quantitatively. Obviously, they can decrease 
the number of flipons and limit their motion. 



B. Collective spin dynamics 

From Eqs.CHll-ITnjl in the Fes crystal at i?^ = 3.8 T 
(Ao « 0.35 K) one finds Tm - 0.1 if for the field com- 
pensated sample {Tm ~ 0.11 K with no field compen- 
sation). The corresponding temperature for the Mni2 
crystal at = 7.0 T (A^ « 0.36 K) is Tm 0.13 K. 
At T < Tm the average distance between fiipons is large, 
and so the collective multi-pair flip-fiop processes are es- 
sentially "frozen" and correlated clusters of fiipons can 
not appear. 

At Tm < T < /S.o/kB the collective multi-pairs pro- 
cesses are unfrozen and, when the temperature increases, 
the whole hierarchy of correlated clusters of increas- 
ing size n (the number of involved fiipons, max{n} = 
Nex(T)) can, in principle, appear. The motion of flipons 
leads to the suppression of correlations. The spectral 
diffusion effect in the limit Aq >> {Wd,Eo} is rather 
weak (Section llll C|) and the spectral diffusion dephas- 
ing time tJ is longer than the motional dephasing time 
T™ . If at certain temperature the cluster dephasing time 
T™''(T) = min{T^"(T),Tj(r)} is shorter than the cluster 
correlation time tdT), creation of the correlated clus- 
ters at average distance Rex{T) is virtually impossible. 
Instead, only short-living (with t^" < tc < Tnd) corre- 
lations between molecules at distances R < Rex can be 
present. 

If fiipons can move, the correlated clusters at aver- 
age distance Rex can appear only at temperatures T ^ 
Tc > Tm (Eqs.lEI 123)7 when r™ becomes longer than 
tc and Rcx{T) becomes shorter then Rc fEa. (|24|l '). In 
Fes (for H"" = 3.8 T) and in Mni2 (for i?^ = 7.0 T) 
crystals this may happen already at T 0.25 K (as- 
suming p — q ~ s — 1/3 and p = 2/3, see Section Fill Cll . 
Note, however, that in the Feg crystal for example, where 
all lattice constants are different, only fiipons oriented 
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(propagating) along the same crystallographic axis can 
create correlated cluster. At the same time, flipons have 
larger probability to propagate along the axis with the 
shortest inter-molecular distance (i.e., with larger AJJ). 
In this case a quasi-ld motion of flipons is more proba- 
ble than a 3c? one. Nevertheless, if at i > flipons will 
start to change their orientation, this process will speed 
up decorrelation. 

If the cluster dephasing time r™'' is tc « r„p 
(Eq|2Jl, during the time interval t = r„p the correlated 
clusters can be created and destroyed (roughly) ~ Tnp/tc 
times. During the cluster life-time r™* all molecules, 
belonging to cluster, can make ^ A^Jtc/^ coherent os- 
cillations (|C7|) . At T ^ 0.25 K the correlated clusters 
can "reappear" ~ 50 times in Fes at = 3.8 T (for 
a cluster of flipons oriented along the same axis) and 
- 10 times in Mni2 at 77^ = 7.0 T. At these fields 
one may expect ~ 30 oscillations (|C7|) in both Feg and 
Mni2. All these estimations do not take into account 
the coherence optimization strategySS and we again as- 
sume p = q = s = 1/3 and p = 2/3. Note also that the 
phonon decoherence rate increases with temperature 
(see EqQJ, but at ksT < Ag this increase is slow. 

In the case of Mni2, in the NPC-window for major 
species the spin dynamics of molecules belonging to mi- 
nor species is incoherent already. Due to the difference 
in Ai, the molecules of major and minor species cannot 
create correlated flip-flop pairs between each other. This 
results (i) in a decrease in the number of flipons; (ii) in 
partial localization of flipons (depending on the fraction 
of impurities) ; and (iii) in randomization of processes due 
to interactions between molecules belonging to different 
species. The last effect gives rise to the incoherent pair 
processes, leading to suppression of coherence. 

The larger the concentration of impurities with inco- 
herent internal dynamics, the smaller the probability for 
correlated clusters to appear. However, at T > Tc the 
sample can become covered by correlated clusters of the 
sizes smaller than the average distance between impuri- 
ties Rim- For example, in Mni2 at = 7.0 T and 
T = 0.25 K one gets Rex ~ 3a. If Rim }t 10a, at these 
values of field and temperature the correlated clusters of 
the radius 3Rex can, in principle, appear. 

If AJj < S^Jj for most of the nearest-neighbor 
molecules (i.e., no flipons), the correlated clusters com- 
posed of lengthy flip-flop pairs with Aff{R ~ Rex) << 
Ajy can appear already at Tm < T < Tc- The asym- 

(2) 

metry for two resonant pairs in such a cluster is 

(2) 

Aff ^ Ajf and the cluster life-time is ~ tc- This 
means that all resonant pairs will be able to make only 
one corelated flip-flop transition before correlations will 
be suppressed (if tc « T„p, such a cluster can reap- 
pear Tnp/tc times). If Ayj > , but for some reason 
/9 ^ (s ^ 1), the correlated clusters of immobile flipons 
can appear. Because of the weakening of the spectral dif- 
fusion effect in the limit Aq >> {Wd^Eq], the dephas- 
ing time for these clusters can be limited only by r„p and 



the number of oscillations ljC7|) can be limited only by 
Q^A^f/Ao- 



V. SUMMARY 

In the present work the internal dynamics of 
a temperature-equilibrated crystalline sample of the 
dipole-dipole interacting molecules with the central spins 
Si has been studied in the coherence window for nuclear 
spin and phonon degrees of freedom. 

At large external transverse magnetic fields the tun- 
neling matrix element Aq (Section ^ and Appendix 
increases, whereas both the half-width of the dipolar bias 
distribution Wd and the half-width of the hyperfine bias 
distribution Eo decrease (Section [HI and Appendix IB)). 

At a certain value of the transverse field, the coher- 
ence window for phonon and nuclear spin-mediated de- 
coherence (the NPC-window) opens up (Sections Q and 
Inj) . Outside of the NPC-window the spin dynamics is 
incoherent. If in the whole NPC-window the average 
strength of the dipole-dipole interactions between the 
nearest-neighbor molecules Vdd or Wd are larger than 
Ao, the spin dynamics is also incoherent. In the opposite 
limit the coherent spin dynamics is possible. 

In the limit Ao > Vdd and if the effective matrix 
element A^™, describing transitions between two fiip- 
flop states of the nearest- neighbor molecules, Eq.®, 
is large compared to the asymmetry of these two 
states, Ea. H10|l . the spin correlations between the nearest- 
neighbor molecules lead to the creation of resonant fiip- 
flop pairs fSection llll Al and Appendix 1^ . Such resonant 
pair experiences oscillations ljC7|l between states | fl-JJ-) 
and I JJ-fl-) of two molecules with frequency ^ ^z/- 

If AJJ > £^fj for the most pairs of the nearest-neighbor 
molecules, the resonant flip-fiop transitions can "prop- 
agate" in the crystal, involving more and more new 
molecules (one molecule in a pair remains in its ground 
state but another nearest-neighbor molecule creates new 
resonance with the excited molecule). This "mobile" 
magnon-like process (a spin-excitation) between the 
states of two involved nearest-neighbor molecules in our 
work is called "fiipon" (Section IIII A|l . The number of 
flipons is limited by the number of excited molecules 
NcxiT)- 

At T < Tm = max{Tiei,rp?,}, Eas. fTM^ . the dis- 
tances between flipons are long and the correlations be- 
tween them are unimportant. At T > Tm the correla- 
tions between flipons become crucial and at certain con- 
ditions can lead to the creation of correlated clusters of 
flipons (Section IIIIB|I . Each cluster represents a corre- 
lated group of molecules experiencing coherent oscilla- 
tions ljC7|l between their lowest states | ff) and | ||.). 

The flipons motion and the spectral diffusion process 
result in the suppression of correlations (Section IIII C|l . 
If AJj >> for most pairs of the nearest-neighbor 
molecules (Section IIV Ap and flipons can move, the cor- 
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related cluster can appear only at T ^Tc, Ea. (|25|l . The 
average inter-flipon distance in this case is given by the 
solution of Ea. H24|l and the cluster dephasing time (i.e., 
its life-time) r^^ is ~ tc (Eq.^). Molecules that be- 
long to the cluster can make ^ AJ^tc/fi oscillations (|C7|) 
before the correlations will be suppressed. During the 
total phonon/nuclear spin coherence time r„p, Eq.(|2Jl, 
the coherent clusters can "reappear" ^ Tnp/tc times. At 
T < Tc only random short-living {t < tc) correlations 
within small groups of flipons, separated by distances 
R < RexiT), Ea. ((TT|l . can appear. 

The smaller the effective flipon diffusion coefficient D / , 
Ea. (jD4p . the longer the cluster life-time t^"- If, for 
some reason, the flipons are localized {Df 0), the 
correlated clusters of immobile flipons can appear even 
at T < Tc- Their life-time can be limited only by T„p, 
and the number of oscillations (|C7|) can be limited only 
by Q^A]]/Ao. If A]] < the correlated clusters 

of "lengthy" resonant flip-flop pairs (molecules in these 
pairs are separated by the distance R ^ R^x ) can appear 
also at T < Tc. Their life-time is limited by tc- On aver- 
age, all pairs in these clusters will be able to make only 
one flip-flop transition before correlations will be sup- 
pressed. If tc < Tnp, these clusters can reappear several 
times. 

It is worth mentioning also that various systems allow 
the coherence optimization strategy (orientation of exter- 
nal transverse field in a plane, chemical replacement of 
isotopes, etci2£) to be applied to get longer spin/phonon 
coherence time-interval Tnp or to shift the coherence win- 
dow down to lower values of the transverse field. 

This concludes our study of the collective spin dynam- 
ics of a temperature equilibrated sample in the coherence 
window for phonon and nuclear spin-mediated decoher- 
ence. The presented analysis can be applied to any crys- 
talline nanomagnetic insulator composed of the central 
spin S molecules, and can be useful in the search for 
magnetic systems showing the spin coherence and collec- 
tive phenomena. The analysis for the induced dynamics 
(when system at t = is prepared in some initial state) 
will be presented separately. 
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This work is supported by Russian grant RFBR 04-02- 
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APPENDIX A: TWO-LEVEL SYSTEM 

The effective Hamiltonin of a biased two-level system 
(TLS) has the form: 

ifTL5 = -AoT^ - (Al) 



where f^,f^ are the Pauli matrixes multiplied by 2; Aq 
is the tunneling matrix element; and ^ is the asymmetry 
between two states (i.e., the longitudinal bias). One can 
easily solve this problem for eigenfunctions: 

= u\^)+v\i)- \i^)^~v\ ^)+u\ i); 
{u,v) = ((e±e)/2e)i/2; e = ^FTAf , (A2) 

where the corresponding energies in states | fl"), | -IJ.) are 
given by -Ef.^ — and 

IT)=Q; U)=(J). (A3) 

If at time t = system was in state | t), the probabilities 
to find system at time t in states | "f) or | |) are 

P„ = l-^sin2(et/n); P^^ = ^ sm\et/h). (A4) 

This describes the oscillations with frequency e between 
states I t) and | |). In the limit Aq << ^ the oscillations 
are suppressed since their amplitude is A^/^^ << 1. 

APPENDIX B: METHOD OF CALCULATIONS 

(1) Wd{H-^). To obtain the transverse field be- 
havior of the dipolar bias distribution half-width Wd in 
the crystals of Fe^ and Mni2 molecules, two clusters 
of different crystal symmetry are used, (a) The Fes 
crystal. The cluster for the Fes system contains 50'^ 
unit cells arranged in a triclinic lattice array with lat- 
tice parametersi^O a = 10.522(7) A; 6 14.05(1) A; c = 
15.(1) A with angles a = 89.90(6)°; /3 = 109.65(5)°; 7 = 
109.27(6)°. Each unit cell of volume K ~ 1969 A^ 
contains eight spin-5/2 Fe~^^ ions, correctly positioned 
and orientedj^fl*^ (b) The Mni2 crystal. The cluster 
for the Mni2 system contains 50'^ unit cells, arranged 
in a tetragonal lattice array with lattice parameters*^ 
a ^ b ^ 17.1627(6) A; c = 12.2880(4) A with angles 
a = /5 = 7 « 90°. Each unit cell of volume K « 3619.5 
A'^ contains twelve spins: four spin-3/2 Mn''+ ions in the 
inner shell and eight spin-2 Mn^"*" ions in the outer shell, 
correctly positioned and oriented.— 

The distributions of the dipolar bias fields and energies 
in the cluster are calculated taking into account all in- 
ternal spins s^f^ of each molecule (Si = J2p p = 8 in 
Fes a-iid p = 12 in Mni2). The internal molecular spins 
sf'^ cannot flip independently - each molecule changes 
its total spin orientation as a rigid object. Initially, all 
molecules in the sample are oriented along the easy axis, 
either at random with projections Sf = zLSi (for depolar- 
ized sample with initial magnetization M = 0), or with 
projections = +Si (for polarized sample with M = 1). 

To obtain the average longitudinal bias field acting 
on i-th molecule, we calculate the longitudinal fields 
h'l'"^\ created by all internal spins of all molecules in the 
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sample at each p-th internal spin s^^-* of i-th molecule. 
Then, the average longitudinal field at j-th molecule is 
Hi = 'Epis'f^/Si)hf^K The calculation of the average 
transverse field is similar. 

First, we calculate the longitudinal and transverse 
fields in the sample, assuming S^ = Si and Sf = 
for each molecule. Knowing internal and external lon- 
gitudinal and transverse fields at each molecule from the 
first step, we calculate A^, S'f and Sf by means of ex- 
act diagonalization of the molecular Hamiltonians H5I6|I . 
Obtaining Sf and Sf for each molecule, we repeat calcu- 
lation of fields in the sample. Using these new fields, we 
recalculate A^, Sf and Sf and so on. This iteration pro- 
cedure converges and, depending on the value of external 
applied field, it is enough to make 10 — 20 iterations to 
obtain a final result. 
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FIG. 5: Total distributions of the internal dipolar fields in tfie 
crystal of Fes molecules (cluster of 50'' molecules arranged in 
a triclinia lattice). The curves are: solid line - the transverse 
dipolar fields distribution for M = 1 (polarized sample) at 
Hfxt = 0; inverse triangles - the longitudinal dipolar fields 
distribution for = 1 at H^^t — 4.8 T; squares - the longi- 
tudinal dipolar fields distribution for M = at H^^^ = 4.8 T 
(depolarized sample - states | -ff) and | ij.) are equipopulated); 
diamonds - the same as squares, but for the model depolarized 
sample (all molecules are in states | -(]■), but Sf /\Sf\ — 0); 
circles - the longitudinal dipolar fields distribution for M = 1 
at H^^t = (polarized sample); triangles - the longitudinal 
dipolar fields distribution for M = at H^^t ~ (depolarized 
sample) . 

To show that this iteration procedure converges, we 
present Fig[Sl where the total distributions of longitudi- 
nal and transverse fields in the crystal of Fes molecules 
are plotted. In a completely polarized sample (M = 1) 
at zero external transverse field Hf^^ the " x" component 



S ZijXij 



/r% (5^ = and 



of the dipolar field is Hfi 
S' = S). At large Hf^t {S' and S'* 
component of the dipolar field is iJ| — ^ju^j^^j / , 
the z-fields distribution at large Hf^^ and the x-fields 
distribution in a polarized sample {M 



S) the "z" 
SxiiZij/rl,. Thus, 



should be the same. This is what one can see in Fig|S| Af- 
ter the iteration procedure the longitudinal fields distri- 
bution for 7\f = at i?ea;t — 4.8 T almost coincides with 
the transverse fields distribution for M = 1 at Hf.J.^ — 0. 

Note also that the bias field distributions in the model 
depolarized sample (all molecules are in states | fl-), but 
^iSf /\Sf\ = 0; shown by diamonds in Fig [51) almost 
coincide with the distributions in the depolarized sam- 
ple (states I i\) and | JJ-) are equipopulated; shown by 
squares). Moreover, the bias distributions in a polar- 
ized (shown by inverse triangles) and depolarized samples 
(shown by squares and diamonds) at large fields almost 
coincide as well. All this means that in order to un- 
derstand the field dependence of important parameters 
at large transverse fields, it is sufficient to calculate all 
these parameters only in the model depolarized sample. 

Repeating the iteration procedure described above for 
each value of external transverse field , we obtain: 1) 
the distributions of longitudinal and transverse dipolar 
bias fields and energies; 2) (5'^)(i?^) and (5'^)(7J^) (the 
sample average absolute values of z and x projections of 
total spin Si); and 3) WoiH'^)- All these calculations 
can be done for any degree of initial polarization M and 
for any degrees of populations of states | ff) and | -IJ-) ; the 
external transverse field can be applied in any direction 
in the x — y plane. Finally, we would like to mention that 
the results obtained for cluster of 40^ molecules do not 
change with further cluster size increase. 

(2) Eo{H-^). All (longitudinal) hyperfine couplings 
{iJl } >tn^t34 i between the Fe (Aln) electronic spins s^^-* 
(5 = X]p^^') the nuclear spins {Ik} are assumed 
dipolar (with the exception of the nuclear spin of any 
Fe^^ isotope in the Fes molecule and of the nuclear 
spin of any MiT'^ nucleus in the Mni2 molecule). The 
strength of the contact hyperfine interaction between the 
Fe electronic spin and the Fe^"^ nuclear spin is known4^ 
(the nuclear spin of the Fe^^ isotope is zero; the standard 
Fes molecule contains > 97% of the Fe^^ isotope). The 
strength of the hyperfine interaction between the Mn^^ , 
Mn^+ electronic spins and the Mn^^ nuclear spin can be 
extracted from the recent NMR measurementst^ 

Knowing the transverse field dependence of (S^) and 
(S^) and the positions and moments of all nuclear spins 
and i^e^+ {Mn"^^, Mn^+) ions in a molecule,— one 

can calculate all the couplings {uj^} and the half- width 



Eo{H^) = iZLiih + l)/fe(c^i')V3)'/' of the hyperfine 
bias energies distribution. All the necessary details of 
calculations can be found in htcraturo^'^i-^^i^^i^-'^i^^i^^). 



1) at Hf^t = 



APPENDIX C: TWO COUPLED TLS 

Consider an ensemble of the interacting two-level sys- 
tems and choose any two coupled systems described by 
the Hamiltonian: 

H = h'.+h'. + J^V^^t^t^; h[ = -A,ff-^:ff . (CI) 

a, 13 
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Here both do not include the bias arising from inter- 
actions with the second TLS (the bias is created by aU 
other TLS in the ensemble), and the last term describes 
the interactions between two systems. It is convenient 
to rewrite the Hamiltonian (|C1|) equivalently, adding the 
contribution coming from the second involved TLS to 
each ^j. In what follows, we limit our consideration only 
by {a,/3} = {x, z}. The term describing the bias created 
by the second TLS at the first one is y^^ff f| + y^^ff 
and the resulting Hamiltonian becomes: 



H = H,+H2 



V 



XX ;:lx ±.x 



(C2) 



where in both Hi and H2 the asymmetries now contain 
the contributions from all the TLS in the ensemble. 

Calculating matrix elements of the Hamiltonian H in 
the representation HA2|I . one gets (the states are ordered 
as{|5i52)} = {|M),l W,l 



H=Ha 



Vdd] 



Ho 



/ -£i - £2 

-£1+62 

£1 - £2 

V £i+£2 



(C3) 



(C4) 



where Ei — \J Q + and 



with A, B, C, D given by 



A 


C 


D 


B 


C 


-A 


B 


-D 


D 


B 


-A 


-C 


B 


-D 


-C 


A 



(C5) 



A = 


£l£2 


^..AiA2 

£l£2 




B = 


Ai Ap 

£1^2 


_|_ yXX 

£ie2 




C = 


yZzCl^2 ^ 

£l£2 


£l£2 




D = 


£l£2 


yXX Cl^2 

£l£2 


(C6) 



In the limit A^ >> matrix elements C and D are small 
in comparison with A and B. 

Note that states | f^) and | JJ.) are the eigenstates of 
the Hamiltonian HA1|I and in both TLS only transitions 
between these states (but not between states | t) and | |) 
(|A3|) 'l are considered. Two central states | ff-IJ.} and | JJ-fl-) 
(from now on we call them the " flip-flop states" ) are sep- 
arated from the two remaining states | ff-ll) and | JJ.JJ.) by 
the energy gaps > A^ and in the limit A^ > > { |fi | , | V"'^ \ } 
the effect of these two remaining states on the flip-flop 
transitions | f|~-|J.) <^ | JJ-Tt) is small (within a second-order 
perturbation theory the corrections to the flip-flop ma- 
trix elements are ~ (F"^)^^?/Af). Therefore, in this 
limit two flip-flop states of a pair can be considered as 



an effective two-level system. The coefficient A in (jC6|) 
plays a role of the mean energy for this effective TLS. 

In the limit A^ >> the tunneling matrix 

element describing the flip-flop transitions is given by 
\B\; the energy change during these transitions is 2|£i — 
£2! << Ai. Then the probability to flnd system at time 
t in state | Jlfl') if at t = it was in state | fl-JJ.) is: 



P< 



m)m) 



A^ 



(C7) 



where the frequency of oscillations, the tunneling matrix 
element connecting states | ff-IJ.) and | JJ-fT) and the asym- 
metry between these two states are given by: 



E 



ff 



ff ■ 



A2 • 



\B\; ^ff ~\£i-£2\ 



(C8) 



(C9) 



If^// < A//, the amplitude of oscillations between states 
I ffjj-) and I JJ-lt) is A^^/i?|^ ~ 1. Since aU other transi- 
tions are described by the matrix elements ^ \V"^\£,i/ Ai 
and the energy change during these transitions is > A^, 
in the limit A^ >> IV"^]} all these transitions are 
not in resonance and their amplitude is small. 

One note is in order here. Applying this model at low 
temperatures to the central molecular spins Si described 
by the Hamiltonians of magnetic anisotropy (|5I6() . one 
should take into account that the spin projection 
along the transverse magnetic flcld has the same sign 
in both exacts lowest states | ff) and | JJ.) of jSJ) and 
The eigenstates of these Hamiltonians can be easily ob- 
tained by the exact diagonalization method. 



APPENDIX D: FLIPON "RANDOM WALKS" 

Consider quasi-ld (i.e., along one crystallographic 
axis) motion of one flipon (the flipon center of mass 
moves in a "fictitious" lattice whose sites are placed di- 
rectly between the nearest-neighbor sites of a real lattice) 
and let us approximate this motion by a " random walks" 
model. In this approximation we assume that flipon can: 
(i) make a jump to the right with the probability p; (ii) 
make a jump to the left with the probability q; and (iii) 
stay at the same site with the probability s = 1 — p — q. 
The probability to find a " walk" with the Kr lattice steps 
to the right and Ki steps to the left from total N steps is 
given by the polynom distribution 



PN{Kr,Ki) 



N\ pKrqKigN-K,.-Ki 
Kr\Kl\{N ~ Kr - Kl)V 



(Dl) 



E PNiKr,Ki) = {p + q + s)^ ^1. (D2) 



Kr=0 K,=0 
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The corresponding displacement for such walk is r = 
Kj. — Ki. The average displacement (?) and dispersion 
{{Ar)^) = (r^) — (r)^ can be easily calculated 

{r)^N{p~q)] {{A¥)^) ^ 4:Npq + Ns{l- s) (D3) 

and for p — q — s~l/3 one gets (r) = and 
((Ar)^) = 2N/3. In our case the time t = NAt is mea- 
sured in units of At = fi/AJ'^(d) and the distance (the 
displacement) r = rAr is in units of Ar = a (these pa- 
rameters describe the flipon elementary "jump"). The 
distribution Ppf{Kr, Ki) defines the normalized distribu- 



tion PAr(r) {r = Kr - Ki G [-7V, TV]). For large values 
of N , this distribution transforms into the Gaussian one 
Pi{r,t) = exp(-rV4£)/t)/(47rD/i)i/2, where Df is the 
flipon effective diffusion coefficient 

Df = pa^A]]/2n; p = 4pq + s{l - s) (D4) 

and {r^{t)} = 2Dft. The Gaussian Pi(r, i) gives the 
probability to find a flipon at time t at the distance r 
from its t = position. 
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